Phase Transitions of the Flux Line Lattice in High- Temperature 
Superconductors with Weak Columnar and Point Disorder 
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We study the effects of wealc columnar and point disorder on the vortex-lattice phase transitions 
in high temperature superconductors. The combined effect of thermal fluctuations and of quenched 
disorder is investigated using a simplified cage model. For columnar disorder the problem maps 
into a quantum particle in a harmonic -I- random potential. We use the variational approximation 
to show that columnar and point disorder have opposite effect on the position of the melting line 
as observed experimentally. Replica symmetry breaking plays a role at the transition into a vortex 
glass at low temperatures. 
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There is a lot of interest in the physics of high tempera- 
ture superconductors due to their potential technological 
applications. In particular these materials are of type II 
and allow for partial magnetic flux penetration. Pinning 
of the magnetic flux lines (FL) by many types of disor- 
der is essential to eliminate dissipative losses associated 
with flux motion. In clean materials below the supercon- 
ducting temperature there exist a 'solid ' phase where 
the vortex lines form a triangular Abrikosov lattice [Q. 
This solid can melt due to thermal fluctuations and the 
effect of impurities. In particular known observed tran- 
sitions are into a flux liquid at higher temperatures via 
a melting line (ML) |2], and into a vortex glass at low 
temperature 1^, jj], [5| in the presence of disorder- the 
so called entanglement line (EL). |l[] 

Recently the effect of point and columnar disorder on 
the position of the melting transition has been measured 
experimentally in the high- Tc material Bi2Sr2CaCu20% 
Q . Point disorder has been induced by electron irradia- 
tion (with 2.5 MeV electrons), whereas columnar disorder 
has been induced by heavy ion irradiation (1 GeV Xe or 
0.9 GeV Pb). It turns out that the flux melting tran- 
sition persists in the presence of either type of disorder, 
but its position shifts depending on the disorder type and 
strength. 

A significant difference has been observed between the 
effects of columnar and point disorder on the location of 
the ML. Weak columnar defects stabilize the solid phase 
with respect to the vortex liquid phase and shift the tran- 
sition to higher fields, whereas point-like disorder desta- 
bilizes the vortex lattice and shifts the melting transition 
to lower fields. In this paper we attempt to provide an 
explanation to this observation. The case of point de- 
fects has been addressed in a recent paper by Ertas and 
Nelson § using the cage-model approach which replaces 
the effect of vortex-vortex interactions by an harmonic 
potential felt by a single vortex. For columnar disorder 



the parabolic cage model was introduced by Nelson and 
Vinokur [8] . Here we use a different approach to analyze 
the cage-model Hamiltonian vis. the replica method to- 
gether with the variational approximation. In the case of 
columnar defects our approach relies on our recent anal- 
ysis of a quantum particle in a random potential |^ . We 
compare the effect of the two types of disorder with each 
other and with results of recent experiments. 

Assume that the average magnetic field is aligned along 
the 2;-axis. Following EN we describe the Hamiltonian of 
a single FL whose position is given by a two-component 
vector r(z) (overhangs are neglected) by: 



\2\dz 



\/(z,r) + ^r2 



(1) 



Here e — eq/t^ is the line tension of the FL, 7^ = 
Tnz/mx^ is the mass anisotropy, eo = (<f>o/47rA)^, (A be- 
ing the penetration length), and /i « £0/^0 the effective 
spring constant (setting the cage size) due to interactions 
with neighboring FLs, which are at a typical distance of 
oo — \/^n/B apart. 

For the case of columnar (or correlated) disorder, 
V{z, r) — V{r) is independent of z, and 



{V{r)V{r')) 
where 

r(2) 



-2/((r 



cr/2)=geiesP{r-r'), (2) 



(r - r') « 1/(27:^) exp(-(r - r')' /2e) 



(3) 



and ^ is the vortex core diameter. The dimensionless pa- 
rameter g is a measure of the strength of the disorder. 
On the other hand for point-disorder, V depends on z 
and§ 



{Viz, r)Viz', r')) = Ae^^^sf (r - r')5{z - z') 



(4) 



The quantity that measures the transverse excursion 
of the FL is 



1 



ul{i)^{\v{z)-v{z + (:)\^) A 



(5) 



Let us now review the connection between a quantum 
particle in a random potential and the behavior of a FL 
in a superconductor. The partition function of the for- 
mer is just like the partition sum of the FL, provided one 
make the identification IHl 



pn^ L, 



(6) 



Where T is the temperature of the superconductor and 
L is the system size in the z-direction. f3 is the inverse 
temperature of the quantum particle. We are interested 
in large fixed L as T is varied, which corresponds to high 
f3 for the quantum particle when h (or alternatively the 
mass of the particle) is varied. The variable z is the so 
called Trotter time. This is the picture we will be using 
for the case of columnar disorder. 

For the case of point-disorder the picture we use is 
that of a directed polymer in the presence of a random 
potential plus an harmonic potential as used by EN. 

The main effect of the harmonic (or cage) potential 
is to cap the transverse excursions of the FL beyond a 
confinement length £* « ao/7. The mean square dis- 
placement of the flux line is given by 



(7) 



The location of the melting line is determined by the 
Lindemann criterion 



(8) 



where cl ~ 0.15 — 0.2 is the phenomenological Linde- 
mann constant. This means that when the transverse 
excursion of a section of length « ^*becomes compara- 
ble to a finite fraction of the interline separation ao, the 
melting of the flux solid occurs. 

We consider first the case of columnar disorder. In 
the absence of disorder it is easily obtained from stan- 
dard quantum mechanics and the correspondence (|^), 
that when L — s- 00, 



u^T) - ^ (1 - exp(~r v/^) = - e-i) 

^ / f 11 \ / ^ fit. 
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from which we find that 
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(9) 
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When we turn on disorder we have to solve the problem 
of a quantum particle in a random quenched potential. 
This problem has been recently solved using the replica 
method and the variational approximation Let us 
review briefly the results of this approach. In this ap- 
proximation we chose the best quadratic Hamiltonian 
parametrized by the matrix Sab{z — z'): 



( dz ( dz'^Sab{z - z')Ya{z) ■ rb{z'). (11) 
Jo Jo „ ■ 



Here the replica index a = 1 . . . n, and n — s- at the end 
of the calculation. This Hamiltonian is determined by 
stationarity of the variational free energy which is given 
by 



In / [dr]exp(-/i„/T), (12) 



where _ff„ is the exact rt-body replicated Hamiltonian. 
The off-diagonal elements of SabCan consistently be taken 
to be independent of z, whereas the diagonal elements 
are z-dependent. It is more convenient to work in fre- 
quency space, where lo is the frequency conjugate to z. 
LUj — (27r/L)j,with j = 0, ±1, ±2, .. ..Assuming replica 
symmetry, which is valid only for part of the tempera- 
ture range, we can denote the off-diagonal elements of 
sab{i^) = (1/T) J^dz e'^^ Sa6(z), by ^(c^) = ^5^,0- De- 
noting the diagonal elements by Sd(w), the variational 
equations become: 



L^,\2T 2T ^ 1 



(13) 



Sd{uj) = s 



T 



dC (1 - e*"^) X 



(14) 



here / '{y) denotes the derivative of the "dressed" func- 
tion f{y) which is obtained in the variational scheme from 
the random potential's correlation function f{y) (see eq. 
(0)), and in 2+1 dimensions is given by: 



fiy) 



2c2 



94^ 



4tt ^2 



y 



(15) 



The full equations, taking into account the possibility of 
replica-symmetry breaking are given in ref . . In terms 
of the variational parameters the function Uq(£*) is given 

by 



2T ^ 1 - cos{uj'£*) 



(16) 



This quantity has not been calculated in ref. There 
we calculated (^r^{0)') which does not measure correla- 
tions along the z-direction. 
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In the limit L ^ oo we were able to solve the equations 
analytically to leading order in g. In that limit eq. jl^ ) 
becomes (for ut ^ Q) : 



SdH = -/ "(&o) 



T 



d(;{l — cos(a;<;)) 



x(/'(Co(0)-/'(M), 



with 



Co{^) = 2T 



du! 1 — cos(cl'<;) 



(17) 



(18) 



and bo given by a similar expression with the cosine term 
missing in the numerator of eq. (|l8|). 
Defining 



(19) 
(20) 

(21) 



T = T / ^/7Ji, a = T /(^^ + r), 
/i(a) = l/(l-a)-(l/a) log(l-a), 

^ OO 

/2(a) = -^(fc + l)«Vfc' 

k=l 

a" = /i(«)//2(a), A = -J"{t) fUa)/Ma)/ii, (22) 
Soo = /"(t) (4 + /i(a))^, (23) 

a good representation of Sd{u!), {oj ^ 0) with the correct 
behavior at low and high frequencies is 



Sdiuj) = Soo + An/ (e uj"^ + a'^n). 



(notice that this function is negative for all w). Substi- 
tuting in eq. (|l^) and expanding the denominator to 
leading order in the strength of the disorder, we get : 



,l{e) = Co(V^) - r(l - A/{a' if/^i) 
x(l-e-^/^')+Tyl/(a(a2-l)V) x 



(1 



)+t/(2^) X {soo+A/ia^-1)) 



X (1 _ e-^/^* - (f/r) e''/''). 



(25) 



In order to plot the results we measure all distances in 
units of ^ , we measure the temperature in units of eq^j 
and the magnetic field in units of $o/?^ ■ We observe 
that the spring constant fi is given in the rescaled units 
by B and oq = 1/VB. We further use 7 = 1 for the 
plots. 

Fig. 1 shows a plot of \/u^j£*)/ao vs. T for zero dis- 
order (curve a) as well as for g/2-K = 0.02 (curve b). We 
have chosen B = 1/900. We see that the disorder tends 
to align the flux lines along the columnar defects , hence 
decreasing u^iT) .Technically this happens since ^d{'^) is 
negative. The horizontal line represents a possible Lin- 
demann constant of 0.15. 

In Fig. 2 we show the modified melting line Bm{T) in 
the presence of columnar disorder. This is obtained from 
eq. (||) with cl = 0.15. We see that it shifts towards 
higher magnetic fields. 



For T < Tc « (eoC/7)[3^eo/(167^VC^)]^/^ there is a so- 
lution with RSB but we will not pursue it further in this 
paper. This temperature is at the bottom of the range 
plotted in the figures for columnar disorder. We will pur- 
sue the RSB solution only for the case of point disorder, 
see below. The expression ( p5| ) becomes negative for very 
low temperature. This is an artifact of the truncation of 
the expansion in the strength of the disorder. 

For the case of point defects the problem is equiva- 
lent to a directed polymer in a combination of a random 
potential and a fixed harmonic potential. This problem 
has been investigated by MP who were mainly con- 
cerned with the limit of ^ 0. In this case the varia- 
tional quadratic Hamiltonian is parametrized by: 

^ rL ^ 

1 



(26) 



with the elements of Sab all constants as opposed to the 
case of columnar disorder. 

The replica symmetric solution to the variational equa- 
tions is simply given by : 



(24) ul{£) = 2T 



1 



Sd 



duj 1 — cos{u!£) 

2ti ELo'^ + ii \ ' e'uj^ + fi 



and hence 



ul{£) 



-_r{l-e-'/'')+T Sd I (2m) 



(27) 
(28) 

(29) 



In eq. (|2^) / is the same function as defined in eq. ( p^ ) 
with g replaced by A. As opposed the case of colum- 
nar disorder, in this case Sd is positive and independent 
of bj, and hence the mean square displacement uf^{£*) is 
bigger than its value for zero disorder. Fig. 1 curve c 
shows a plot of ^/ul{l*)/aQ vs. T for A/27r = 0.8. Again 
B = 1/900. For T <TcpK (eoC/7)(7 A/27r)i/3 it is nec- 
essary to break replica symmetry as shown by MP [p^ . 
This means that the off-diagonal elements of the varia- 
tional matrix Sab are not all equal to each other. MP 
worked out the solution in the limit of — s- 0, but it 
is not difficult to extend it to any value of ji. We have 
worked out the first stage RSB solution which is all is 
required for a random potential with short ranged corre- 
lations. The analytical expression is not shown here for 
lack of space. The solution is represented by curve d in 
Fig. 1 which consists of upward triangles. 

The modified melting line in the presence of disorder 
is indicated by the curve c in Fig. 2 for T > Tcp. For 
T < Tcp the so called entanglement line is represented 
by curve d of filled squares. The value of the magnetic 
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field Bm{Tcp) « ($o/C2)(7A/27r)-2/3ci gives a reason- 
able agreement with the experiments. 

The analytical expressions given in eqs. (p5|), (29), 
though quite simple, seem to capture the essential feature 
required to reproduce the position of the melting line. 
The qualitative agreement with experimental results is 
remarkable, especially the opposite effects of columnar 
and point disorder on the position of the melting line. 
The 'as grown' experimental results are corresponding 
to very small amount of point disorder, and thus close to 
the line of no disorder in the figures. At low temperature, 
the entanglement transition is associated in our formal- 
ism with RSB, and is a sort of a spin- glass transition in 
the sense that many minima of the random potential and 
hence free energy, compete with each other. In this paper 
we worked out the one-step RSB for the case of point dis- 
order. The experiments show that in the case of colmunar 
disorder the transition into the vortex glass seems to be 
absent. This has to be further clarified theoretically. We 
have shown that the cage model together with the varia- 
tional approximation reproduce the main feature of the 
experiments. Effects of many body interaction between 
vortex lines which are not taken into account by the ef- 
fective cage model seem to be of secondary importance. 
Inclusion of such effects within the variational formalism 
remains a task for the future. 

For point disorder, in the limit of infinite cage ( /i ^ 0), 
the variational approximation gives a wandering expo- 
nent of 1/2 for a random potential with short ranged 
correlations jl^], whereas simulations give a value of 5/8 
PI . This discrepancy does not seem of importance with 
respect to the conclusions obtained in this paper. An- 
other point to notice is that columnar disorder is much 
more effective in shifting the position of the melting line 
as compared for point disorder in the range of parame- 
ters considered here. We have used a much weaker value 
of correlated disorder to achieve a similar or even larger 
shift of the melting line than for the case of point disor- 



der. The fact that the random potential does not vary 
along the z-axis enhances its effect on the vortex lines. 
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We thank the Weizmann institute for a Michael Visiting 
Professorship, during which this research has been car- 
ried out. 



[1] 
[2] 

[s: 
[4: 
[s: 

[6: 

[7 

[8 

[9: 

[lo; 

[11 



G. Blatter et al, Rev. Mod. Phys. 66, 1125 (1994). 

E. Zeldov et al, Nature 375, 373 (1995); see also H. Pas- 

toria et al. Phys. Rev. Lett. 72, 2951 (1994) 

M. Feigelman, et al, Phys. Rev. Lett. 63, 2303 (1989); 

A. L Larkin and V. M. Vinokur, ibid. 75, 4666 (1995). 
D. S. Fisher, M. P. A. Fisher and D. A. Huse, Phys. Rev. 
B43, 130 (1990). 

T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72, 
1530 (1994); Phys. Rev. B52, 1242 (1995); see also 
T. Nattermann, Phys. Rev. Lett. 64, 2454 (1990). 

B. Khaykovitch et al., Phys. Rev. Lett. 76, 2555 (1996) 
and preprint (1996). 



D. Ertas and D. R. Nelson, preprint, cond-mat/9607142 
(1996) 

D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); 
D. R. Nelson and V. M. Vinokur, Phys. Rev. B48, 13060 
(1993). 

Y. Y. Goldschmidt, Phys. Rev. E 53, 343 (1996); see also 
Phys. Rev. Lett. 74, 5162 (1995) 

M. Mezard and G. Parisi, J. Phys. I (France)l, 809 (1991) 
T. Halpin-Healy and Y.-C. Zhang, Phys. Rep. 254, 215 
(1995) and references therein. 



Figure Captions: Figl: Transverse fluctuations in the 
cage model for (a) no disorder (b)columnar disorder 
(c)point disorder (d)RSB for point disorder. Fig. 2: 
Melting line for (a) no disorder (b) columnar disorder 
(c)point disorder (d) entanglement line for point disor- 
der. 
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Figure 2 
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